Fidelity of a Rydberg blockade quantum gate from simulated quantum process 

tomography 
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We present a detailed error analysis of a Rydberg blockade mediated controlled-NOT quantum 
gate between two neutral atoms as demonstrated recently in Phys. Rev. Lett. 104, 010503 (2010) 
and Phys. Rev. A 82, 030306 (2010). Numerical solutions of a master equation for the gate 
dynamics, including all known sources of technical error, are shown to be in good agreement with 
experiments. The primary sources of gate error are identified and suggestions given for future 
improvements. We also present numerical simulations of quantum process tomography to find the 
intrinsic fidelity, neglecting technical errors, of a Rydberg blockade controlled phase gate. The gate 
fidelity is characterized using trace overlap and trace distance measures. We show that the trace 
distance is linearly sensitive to errors arising from the finite Rydberg blockade shift and introduce a 
modified pulse sequence which corrects the linear errors. Our analysis shows that the intrinsic gate 
error extracted from simulated quantum process tomography can be under 0.002 for specific states 
of 87 Rb or Cs atoms. The relation between the process fidelity and the gate error probability used 
in calculations of fault tolerance thresholds is discussed. 

PACS numbers: 03.67.-a, 32.80.Qk, 32.80.Ee. 



I. INTRODUCTION 

Arrays of laser cooled neutral atoms in optical lattices 
or far-off-resonance optical traps (FORTs) arc promising 
candidates for quantum computing experiments, due to 
their very long decoherence time (up to several seconds) 
for ground state atoms and strong two-atom interaction 
for highly excited Rydberg state atoms. This strong, 
long-range and controllable interaction leads to the so- 
called Rydberg blockade effect in which only one atom 
in an ensemble can be excited into a Rydberg state if 
the ensemble size is smaller than the Rydberg blockade 
radius. Jaksch et al. [l[ first proposed to use Rydberg 
blockade to implement a fast two-qubit controlled-phase 
gate (Cz), which can be converted into a controlled-NOT 
gate (CNOT) using single qubit rotations @]. Soon after, 
various schemes were proposed for fast quantum gates 
with an ensemble 0-^1, entangled state preparation Q, 
quantum algorithms @, , quan tum simulators , and 
efficient quantum repeaters 

Rydberg blockade, the central ingredient of the above 
schemes, has been demonstrated between two individual 
neutral atoms held in FORTs [HI, [l3| , and was used to 
demonstrate a two-qubit CNOT quantum gate [l4[ and 
to generate entangled Bell states with fidelity of about 
0.58 - 0.75 after atom loss correction (0.46 - 0.48 with- 
out atom loss correction) [l4T - [l6l ]. Using an improved 
apparatus [l7| a fidelity of 0.92 for the CNOT probabil- 
ity truth table and 0.71 for Bell state fidelity after atom 
loss correction (0.71/0.58 for the CNOT truth table/Bell 
state fidelity without atom loss correction) was demon- 
strated. These proof-of-principlc results arc promising 
but are still far from simple estimates of Rydberg gate 
errors at the level of E ~ 10~ 3 predicted in [iT [l8j. 

The simple estimates are based on intrinsic errors as- 
sociated with the atomic physics of the states used for 
Rydberg blockade. The essential intrinsic errors are the 



finite lifetime of Rydberg states and the finite strength of 
the Rydbcrg-Rydberg blockade interaction. In addition 
to intrinsic errors, experiments are sensitive to several 
different sources of technical error: spontaneous emis- 
sion from an intermediate level in a two-photon excita- 
tion scheme, magnetic field fluctuations, pulse area er- 
rors, Doppler effects due to finite atomic temperature, 
etc.. We have previously estimated gate errors [iT], EH 
by treating each error source separately in the small error 
limit and adding them together. This provides a good 
estimate for small errors but it is unreliable for larger 
errors, as in the experiments, and does not provide a 
rigorous fidelity measure for the gate operation. 

In this paper we present a more rigorous treatment by 
including all known error sources in a master equation for 
the density matrix evolution (optical Bloch equations). 
We dynamically track the density matrix evolution dur- 
ing the CNOT pulse sequence and then average the re- 
sults over the computational input states for comparison 
with simple analytical estimates. We also extract a rig- 
orous value for the gate fidelity using simulated quantum 
process tomography. The numerical results are in good 
agreement with analytical error estimates when technical 
errors are neglected, and with experimental results when 
technical errors arc included in a Monte Carlo simula- 
tion. Numerical simulations of quantum process tomog- 
raphy with realistic atomic parameters confirm that it is 
in principle possible to reach quantum process errors of 
2 x 10~ 3 for both Rb and Cs atoms. 

In Sec. [II] we present the experimental setup and pro- 
cedures used to demonstrate the CNOT gate and gen- 
erate entangled Bell states. We also enumerate the var- 
ious sources of technical imperfection in a realistic ex- 
periment. In Sec. IIIII we give analytical estimates of 
the intrinsic gate error and present a new pulse sequence 
which removes the leading linear term in the finite block- 
ade shift error. In Sec. IIVI we present a master equation 
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model which includes both the technical error sources 
from [H] and intrinsic errors from Sec. IIHI In Sec. IIV Al 
we compare numerical Monte Carlo simulations with ex- 
perimental results, and demonstrate good agreement. In 
Sec. IIV Bl we perform simulated quantum process to- 
mography of a two qubit controlled-phase gate account- 
ing only for intrinsic errors. This analysis shows that 
in a well designed experiment where technical errors are 
minimized it should be possible to reach low gate errors, 
below known fault tolerance thresholds [lil [2(|, for both 
Rb and Cs. A summary and discussion is presented in 
Sec. H 



II. EXPERIMENTAL PROCEDURES AND 
TECHNICAL IMPERFECTIONS 
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FIG. 1. (color online) a) Experimental geometry. Two single 
atoms are trapped in two FORTs separated by about 9 /tin. 
b) Relevant levels of S7 Rb. To reach the Rydberg state we 
use two-photon excitation with wavelengths of 780 nm and 
480 nm. 



The experimental apparatus and procedures, as shown 
in Fig. [TJ have been described in detail in our recent pub- 
lications [13, [13] ■ For ease in understanding the subse- 
quent analysis we include a brief description of the proce- 
dures, as well as estimates of various sources of technical 
imperfections. 



A. Experimental procedures 

We use FORTs to localize single 87 Rb atoms, which can 
be individually addressed using tightly focused beams 
that are scanned by acousto-optical modulators. The 
FORT beams, propagating along +z, are formed by fo- 
cusing a laser beam with wavelength of Af = 1064 nm to a 
waist (1/e 2 intensity radius) of Wf = 3.4 ^m. We generate 
a linear array of trap sites using a diffractive element with 
the central site's trap depth of Uo/kB = 4.5 mK and trap 
separation of about 9 /im along the x direction. We use 
two sites, one labeled as the control and the other as the 
target, to perform two-qubit quantum logic operations 
and to generate entangled states. A bias magnetic field 
is applied along z, which defines the quantization axis for 
the optical pumping (B z = 0.15 mT) and lifts the degen- 
eracy of the Zeeman sublevels (B z = 0.37 mT) of Ryd- 
berg states during the gate operation. The relevant levels 
of 87 Rb are shown in Fig. [TJd. We use the 5si/2 hyper- 
fine clock states as our qubits |0) = \F = 1, mp = 0) and 
1) = \F = 2,m F = 0), separated by uj 10 /2it = 6.83 GHz, 
and the Rydberg state \r) = \97d 5 / 2 ,mj = 5/2). 

We perform single qubit rotations between |0) and |1) 
using two-photon stimulated Raman transitions driven 
by focusing a er + polarized 780 nm laser with frequency 
components separated by luiq and detuned by 101 GHz 
to the red of the D2 transition [2l[ . Total typical power 
in the two Raman sidebands is ~ 90 /iW with waist of 
Ux/y^g = 7.7 im\, giving a single qubit Rabi frequency of 
Q,g/2iT = 0.56 MHz with ir pulse times of ~ 900 ns and 
peak-to-peak amplitude of better than 0.98 after correc- 
tion for background atom loss of 10% (l7j . 

For coherent Rydberg excitation between |1) and \r) 



we use a two-photon transition with er + polarized 780 
and 480 nm beams [22j ]. They counter-propagate along 
the trap's axial z direction to minimize Doppler broaden- 
ing of the transition. The 780 nm beam is tuned about 
2.0 GHz to the red of the |1) -> \5p 3/2 ,F' = 3) tran- 
sition with typical beam power of 2.4 /iW and waist of 
^x/y^R — 7.7 /im. The 480 nm beam is tuned about 
2.0 GHz to the blue of the \5p 3 / 2 ,F' = 3) -> |r) tran- 
sition with typical beam power of 13 mW and waist 
of (-L> x /y,B = 4.5 /im. This gives a Rydberg red Rabi 
frequency of Or/27t = 118 MHz, a Rydberg blue Rabi 
frequency of Ob/27t = 39 MHz, and Rydberg Rabi fre- 
quency of Q/27T = 1.15 MHz with ir pulse times of 
~ 440 ns and amplitude of 0.92 after correction for back- 
ground atom loss of 10%. 

In order to perform a two qubit CNOT gate we start by 
loading single atoms from a background vapor magneto- 
optical trap into two FORT sites. The trapped atoms 
have a measured temperature of T a ~ 175 /xK using a 
release and recapture method [23[ . Atom detection is ac- 
complished by simultaneously illuminating both trapping 
sites with near resonant red-detuned molasses light, and 
imaging the fluorescence onto a cooled EMCCD camera. 
Detected photon counts arc integrated for approximately 
20 ms. Comparison of the integrated number of counts in 
a region of interest with predetermined thresholds indi- 
cates the presence or absence of a single atom [l2[ • After 
single atoms are loaded in these sites they are optically 
pumped into |1) with efficiency of about 99% [131 us ~ 
ing 7r polarized light propagating along ~x tuned to the 
|5si/ 2 , F = 2) ->• |5pi/ 2 , F' = 2) Dl transition at 795 nm 
and \5si/2,F = 1) — > |5p3/2,F = 2) D2 transition at 
780 nm. This is followed by ground state 7r pulses to 
either or both of the atoms to generate any of the four 
computational basis states. We then turn off the opti- 
cal trapping potentials for about 4 /is, apply the CNOT 
pulses, and restore the optical traps. Ground state 7r 
pulses are then applied to either or both atoms to se- 
lect one of the four possible output states. Atoms left in 
state |1) are removed from the traps with unbalanced ra- 
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diation pressure (blow away light), and a measurement is 
made to determine if the selected output state is present. 
The selection pulses provide a positive identification of 
all output states and we do not simply assume that a 
low photoelectron signal corresponds to an atom in |1) 
before application of the blow away light [HI, H3] • 

Following the above procedures we have obtained the 
CNOT truth table fidelity of F = jTr[U£ e JJ C NOT] = 
0.92 ± 0.06 with [/ideal and Uqnot the ideal and ex- 
perimentally obtained CNOT truth tables. To measure 
the state preparation fidelity, we use the same sequence 
but without applying the CNOT pulses. The computa- 
tional basis states are prepared with an average fidelity 
ofF = 0.97. 

To create entangled states we use ir/2 pulses on the 
control atom to prepare the input states \ct) = 775 (|0) + 

i|l»|l) and \ct) = ^(|0) + »|1»|0>. Applying the 
CNOT to these states creates two of the Bell states 
\Bi) = ^(|00) + |11)) and \B 2 ) = ^(|01) + |10». In 
order to verify entanglement we measured the coherence 
of the |.Bi) state by parity oscillations flil. and ob- 
tained a Bell state fidelity of 0.58 without any atom loss 
correction (0.71 after correction) (l7| . Comparable num- 
bers for the Bell state fidelity were obtained in a related 
experiment (l5| . 



B. Technical imperfections 
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FIG. 2. (color online) Calculated blockade shift B from the 
theory of Ref. [25| and double excitation probability P2 for 
state \97d 5 /2,rrij — 5/2} as a function of relative position 
P(\zi — Z2 1) for 0.37 mT magnetic field. The relative proba- 
bility distribution P(|zi — Z2I) assumes a trapping potential 
U/kB = 4.5 mK, waist w — 3.4 /im of the 1064 nm trapping 
light and atom temperature T= 175 /J.K. 
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blockade shift due to the atomic position distribution, we 
use the theoretical blockade shift curve as a function of 
relative atomic separation \zi — z 2 \ as shown in Fig. [2] 
The position dependent Rabi frequency and AC Stark 
shifts are 



The main technical errors that affect the CNOT oper- 
ation are spontaneous emission of the intermediate level 
during Rydberg excitation, magnetic noise, Doppler ef- 
fects due to finite atom temperature and Rydberg laser 
power fluctuations. 

The spontaneous emission error of the intermediate 
level during a tt excitation pulse can be estimated by 
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linewidth and detuning from the p level. This error is 
about 0.8% for our current experimental setup. 

Magnetic field fluctuations cause transition shifts, giv- 
ing a Rydberg two-photon detuning Ab = (g r iTi r — 
<7i?Tii)^B_B 2 with g r = 6/5, m r = 5/2, g\ = 1/2, and 
mi = for our implementation in Fig. [Tp). We assume 
that the magnetic field fluctuations are Gaussian dis- 
tributed with a standard deviation of <tb = 2.5 x 10~ 6 T. 
This value was found by measuring the decoherence 
time of the hyperfine qubit at two different bias mag- 
netic field strengths [24j. Doppler broadening at finite 
atom temperature T a , also gives two-photon detuning 



with 



A,, the radiative 
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(kji-kB)v with averaged variances of a 2 x = a u 



for both control and target atoms. 



Other technical errors associated with finite atom tem- 
perature are Rydberg blockade shifts, pulse area fluctu- 
ations, and AC Stark shifts (or two-photon detunings) 
due to the atomic position distribution in the FORTs. 
We assume that the atomic position distribution for both 
control and target atoms is Gaussian with variance of 18 1 



n H (r) = Or(0)- 



n B (r) = n B (o)- 



[(l + zV^,«)(l + ^ fl )] 1/4 



[(l + z 2 /^, B )(l + ^/^)] 1/4 



(1) 



(2) 



and 



A AC (r) = A AC (0) 



4IAJ 



(3) 



Here Or(0), ^b(0), and Aac(0) are the Rydberg red, 
blue and AC Stark shifts at trap center, respectively; the 
Rayleigh lengths are L x / VyR = ^w 2 x/y R /\^ and L x / y>B = 

Power fluctuations of the Rydberg lasers will not only 
affect the Rabi frequencies and J7b, but will also af- 
fect the two photon detuning Aac- We assume that the 
power fluctuations of the red and blue lasers are both 
Gaussian distributed with FWHM of 1% and 2%, respec- 
tively, as measured independently. 
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C. Dephasing errors 

The technical imperfections listed in the previous sec- 
tion show up as errors in the measured CNOT probability 
truth table as well as in the fidelity of the output quan- 
tum states. There are additional technical dephasing er- 
rors that do not significantly affect the CNOT truth table 
but strongly impact the fidelity of Bell state generation. 
As was pointed out in Refs. [TBI, [H[ both magnetic field 
fluctuations and atomic motion lead to dephasing of the 
Rydberg state relative to the ground state during gate 
operation because the motion of Rydberg atoms between 
excitation and deexcitation pulses leads to a stochastic 
phase that degrades the entanglement. In the numerical 
simulations described in Sec. [IV] we do not keep track 
of the position dependent phase of the optical fields in 
the evolution Hamiltonians during each blockade pulse 
sequence(Eqs. (|14b[) and (|15a[) below) . Instead we add 
an extra dephasing term 7 p h to the Liouville operators 
of Eqs. (pb|) and (fl5b1 in Sec. IWl 



7 P h 



V(7b) 2 + (td) 



(4) 



where 7b = |Ab|/7i and 7d = Ad are the dephasing rates 
due to magnetic field fluctuations and Doppler effects, re- 
spectively. In the Monte Carlo simulations presented be- 
low these dephasing rates are sampled from distributions 
that are generated with the position and velocity vari- 
ances described above. Both the magnetic field fluctua- 
tions and atom position variations at finite temperature 
also dephase the qubit states by varying the hyperfine 
splitting between them. We model the qubit dephasing 



7oi = \J (7boi) 2 + (ir) 2 



(5) 



where 7301 = wio(-B* + cr B ) - u w (B z ), is the dephas- 
ing rate due to the second-order Zeeman shift of the 
clock transition by magnetic field fluctuations, B z is a 



static bias field, uj\q(B z 



(gs-gi)nBB z 



the ground hyperfine splitting of the clock states, gs,gi 
are electron and nuclear Lande factors, and or is the 
magnetic field fluctuation; 7^ = 1.03 



5 k B T a 
2U 



IS 

the dephasing rate in 1/ms due to atomic motion in 
the FORT, U is the peak FORT potential, 5 /2tt = 
{Uq/Ub) x (1.5 kHz/mK), is the peak differential light 
shift of the FORT and T a is the atom temperature. 

In addition to the above errors, there are also errors 
associated with state preparation due to imperfect op- 
tical pumping and single qubit rotations. These errors 
are at about the 1% level. There is also about 1% atom 
loss due to background collisions before the CNOT pulse 
sequence. 



III. INTRINSIC ERROR ESTIMATES 

Even if all sources of technical error described in Sec. 
IIIBI arc negligible there will still be intrinsic gate errors 



due to the basic physics of the Rydberg blockade interac- 
tion. Intrinsic errors of a Rydberg blockade CNOT gate 
include the decoherence error due to the finite lifetime r 
of the Rydberg state and rotation errors due to imperfect 
blockade. In the strong blockade limit (f2 <C B <C u>io, 
where B is the Blockade shift), the intrinsic gate error 
E\ averaged over the input states in the computational 
basis ({|00>, |01), |10), |11)}) is [H,|l7| 
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The first term in Eq. (|6|) is the Rydberg decay error 
due to the finite lifetime r of the Rydberg state, and the 
second term is the imperfect blockade error. In the limit 
of W10 3> (B, fi) we can extract a simple expression for 
the optimum Rabi frequency which minimizes the error 
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error of 



Oi o t leads to a minimum averaged gate 
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In our experiments r ~ 300 fxs is the radiative lifetime 
of the 97c? 5 /2 Rydberg level, ujio/2tt — 6.83 GHz, and 
f2/27r = 1.15 MHz. In the experimental geometry shown 
in Fig. QJi) a range of two-atom separations, and hence 
blockade shifts, occur. The blockade shift curve shown in 



Fig. [21 is calculated from the theory of Ref. [25[ using a 
trap separation of a; = 8.7 fim and a bias magnetic field 
of B z = 0.37 mT applied along the z axis. Averaging 
Eq. ([6]) over the probability distribution P(|zi — 22 1), 
which is dependent on the trapped atom temperature of 
175 fiK, gives an expected error of E = 8.5 x 10 -3 . The 
corresponding averaged blockade shift from Eq. (|6|) is 
B/2tt = 5.3 MHz. 

It should be emphasized that the average error E\ dis- 
cussed above ignores errors in the phases of the states 
generated by the CNOT gate and therefore corresponds 
to the measurement of a probability truth table. As dis- 
cussed in Ref. [l|, a Rydberg- mediated Cz gate has a 
phase error for the |11) input state of |11) — > — e^|ll) 
with (j> = 7iT2/2B. Averaging over the four computa- 
tional basis states gives an average intrinsic phase error 
of E$ = 7rf2/8B. Including the phase error, we find an 
average intrinsic gate error of 



E-, 



^+8B- 



(9) 



When SI -c B the last term in Eq. ©, dominates over 
the imperfect blockade term in Eq. ©. Using Eq. © 
we can again extract a simple expression for the optimum 
Rabi frequency which minimizes the error 
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FIG. 3. (color online) Pulse sequence for Rydberg Cz gate 
with two-photon 780 and 480 nm excitation of Rb atoms 
following Pi a) and modified sequence to correct finite block- 
ade phase error by rapidly changing the phase of one of the 
Rydberg beams by (f> b). 



in a phase gate 
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Adding a single qubit Z rotation of —<fi to the control 
atom cancels the extra phases, leaving an ideal Cz oper- 
ation. As we show in detail in Sec. IIVBI using simulated 
process tomography, the standard sequence of Fig. [3^ 
gives linear phase errors that impact the trace distance 
fidelity measure. Using the modified pulse sequence the 
linear phase errors are canceled and the trace distance 
error is reduced. 

We emphasize that perfect correction of the phase er- 
ror using this modified gate sequence assumes that B, 
and therefore the phase 4>, are well defined, and do not 
fluctuate. This is not true for thermally excited atoms 
but will be the case for atoms that are in the motional 
ground state of the confining potentials. If the atoms are 
thermally excited we can still cancel the average value of 
<fi, but errors due to fluctuations about the average would 
remain. In the analysis of the intrinsic error limit we as- 
sume that the atoms are in the motional ground state so 
the phase can be canceled exactly. 



IV. NUMERICAL SIMULATIONS 
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Comparing Eqs. (|5J| and in the experimentally rel- 
evant limit of Bt ^> 1 we see that i?2 min 3> Ex min . This 
result seems to imply that Eq. ([5J provides an overly 
optimistic estimate for the gate error. 

Despite this apparently disappointing result the gate 
error can indeed be made to satisfy the (Br) -2 / 3 scal- 
ing of Eq. ([8]) by a modification of the Cz sequence as 
shown in Fig. |3Jd). We note that the only input state for 
which the target atom becomes Rydberg excited is 1 01) . 
By shifting the phase of the |5p 3 / 2 ) — > \r) laser used for 
Rydberg excitation by an amount <j) during the gate op- 
eration we can add a phase of <fi to this state. This results 



In order to accurately simulate the performance of the 
Rydberg gate we integrate the master equation for the 
two-atom dynamics including all known coherent and in- 
coherent rates. A related analysis was performed pre- 
viously for the Rydberg excitation dynamics of a single 
atom |28[ . For each atom we include five atomic states in 
our numerical calculation, which are labeled in the level 
scheme of Fig. [TJd): qubit |0), qubit |1), and reservoir 
level \g) = |5si/2, mp ^ 0) in the 5si/2 ground state; the 
intermediate level \p) = |5p3/2) and the Rydberg level 
\r). With this set of basis states the two-atom dynamics 
are described by density matrices p c t(t) with dimensions 
25 x 25. We take the initial condition to be a separable 
state Pct(0) = p c (0) ® Pt(0), with c/t for control/target 
atoms. We calculate the time evolution by solving the 
master equation 



d-Pct 
dt 



■[H ct ,p ct ] + C c 



(13) 



with H r 



h + L 



H t + hB 



24 
1 

£ c g) J t + 7 C (g) £ tj J t (7 C ) are 5x5 identity matrices, 
and O24 is a 24 x 24 zero matrix. The Hamiltonian H c 
{H{), after making the rotating- wave approximation, and 
the Liouvillc operators C c (C t ) are given in the basis 
{|0),| 5 ),|l)»,|r)}as 
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For Rb atoms we take into account the spontaneous 
emission from the intermediate level \p) to the ground 
level with a decay rate of j p /2ir = 6.07 MHz and the 
corresponding branching ratio of 0.56 to state \g), 0.32 
to state |1) and 0.12 to state |0) as well as the decay 
from the Rydberg state \r) to the intermediate level \p) 
with rate j r /2Tr = 0.53 kHz. A p is the intermediate level 
detuning, A r = Aac(i") + Ab + Ad is the two photon 
detuning (see Sec. IIIBj) . ujiq/2tt = 6.83 GHz is the hy- 
pcrfine ground state splitting; ji r = 7,. + 7 p h is the total 
dephasing of the Rydberg state relative to |1), 7 P h is the 
Rydberg dephasing rate due to magnetic field fluctua- 
tions and Doppler effects as shown in Eq. @; and 701 
is the dephasing of qubit states due to magnetic field 
fluctuations and atomic motion. 

Note that we do not include driving terms in Eq. (|14b|) 
that couple the reservoir level \g) back to \p) and \r). Do- 
ing so correctly would require adding additional Rydberg 
levels with different values of rrij which would increase the 
numerical burden. Since any population in \g) is already 
fully counted as an error, including additional driving 
terms would only reduce the final errors, and our results 
are reliably upper bounds on the gate error. 



A. Monte Carlo Simulations including technical 
errors 

In our numerical calculation, we consider two traps 
aligned along z with separation of 8.7 ^m as shown in Fig. 
[IJi). The atoms in each trap have temperature T a and 
Gaussian spatial probability distribution with variances 
of a x , <T y , and u z as given in Sec. Ill Bl Rydberg Rabi 
pulses are applied to the control or target atoms with 
Gaussian power fluctuations of FWHM of 1% and 2% for 
Rydberg red and blue lasers, respectively. An atom at 
position r = (x,y, z) with velocity v = (v x , v y , v z ) ex- 
periences a Rydberg excitation pulse with effective Rabi 
frequency and two-photon detuning that depends on po- 
sition and velocity, so il = f2(r) and A r = A,.(r,v) as 
described in Sec. Ill Bl We also use a fit to the blockade 
curve of Fig. [5] to account for variations in the two-atom 
interaction strength. Finally we can monitor the effect of 
various error sources by switching them on and off and 



TABLE I. Experimental parameters used in Ref. [l7| as well 
as in the numerical simulations. 



Experimental parameter 


symbol 


value 


FORT wavelength 


A/ 


1064 nm 


FORT waist 


Wx,y 


3.4 fim 


FORT trap depth 


U /k B 


4.5 mK 


FORT separation 


d 


8.7 /urn 


Atom temperature 




175 ftK 


Rydberg level 


|r) 


97d 5 / 2 


Rydberg state radiative lifetime 


r 


320 fis 


Blockade shift at jj,K 


Bo/2tt 


20 MHz 


Rydberg red power 




2.4 jiW 


Rydberg red waist 


w x/y,R 


7.7 ^m 


Rydberg red detuning 


A p /2n 


-2 GHz 


Rydberg blue power 




13 mW 


Rydberg blue waist 


w x/y,B 


4.5 fim 


Rydberg red Rabi frequency 


fi R /27T 


118 MHz 


Rydberg blue Rabi frequency 




39 MHz 


Rydberg Rabi frequency 


Q/2-k 


1.15 MHz 


Magnetic field fluctuation 


OB 


2.5 jiiT 


Rydberg red power fluctuation 




1% 


Rydberg blue power fluctuation 




2% 



comparing the numerical results. 

Numerical simulation of the CNOT gate demonstrated 
in Ref. [l?} proceeds as follows. We start with the initial 
density matrix for both atoms with population of 0.98 in 
either |1) or |0) and population of 0.02 in \g) to account 
for an optical pumping error of 1% and atom loss of 1% 
before the CNOT pulses. Then we perform Monte Carlo 
simulations of the experiment with the actual experimen- 
tal parameters from Ref. pj} as listed in Table |U Next, 
we solve the time evolution of the master equation f| 13[) 
for the Cz pulse sequence of Fig. with 7r/2 pulses 
on the target atom before and after the Cz to give a 
CNOT operation: {(§ ) t , tt c , t gap , (2ir) t , t gap , tt c , (§) t } as 
in [TtJ and average over the four input states ( 1 00) , 1 01) , 
1 10) and 1 11)) to get the averaged gate errors, where (|) t 
is a ground Raman 7r/2 pulse on the target atom, and 
^gap = 500 ns is the switching time between control and 
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TABLE II. CNOT probability truth table and entanglement 
fidelity results. The upper section of the table gives error bud- 
get values used in the numerical simulations. These values are 
based on laboratory measurements and theoretical estimates 
as described in [13]. The lower section gives measured results 
from [l7l ] and numerical values from solving Eq. (|13[) . 





error 




Parameter used in numerical simulation 


budget 




Optical pumping 


0.02 




Atom loss before CNOT pulses 


0.02 




Spontaneous emission 


0.018 




Rydberg decay error 


0.003 




Blockade error at fiK 


0.0004 




Blockade error at 175 i-iK 


0.006 




Doppler Broadening at 175 /j,K 


0.003 




Laser power fluctuation 


0.0001 




Magnetic field fluctuation 


0.0002 






measured 


numerical 




results 


simulation 


Background loss (two atoms) 


0.19 




CNOT trace loss (1 - Tr[p c t]) 


0.01 




CNOT probability truth table 


raw fidelity 


0.74 


0.75 


background loss corrected 


0.91 


0.93 


background & trace corrected 


0.92 


0.93 


Bell state 


raw fidelity 


0.58 


0.54 


background loss corrected 


0.71 


0.67 


background & trace corrected 


0.71 


0.67 



50 evolutions of the master equation using Monte Carlo 
simulation of all the error sources as mentioned before, 
and obtain the final Bell state fidelity F = 0.54 without 
atom loss correction and fidelity F = 0.67 after atom loss 
correction as shown in Table [TT1 which is consistent with 
our measured entangled fidelity of 0.71 in [l7j . We as- 
sume that the atom loss due to collisions with untrapped 
background atoms is independent of the CNOT pulse se- 
quence (about 4 /is) which is much shorter than the trap 
lifetime (several seconds), so the background loss is sim- 
ply considered as a scaling factor for the final fidelity 
without atom loss. 

In Figure SI we compare the numerical simulation re- 
sults for the CNOT probability truth table with the an- 
alytical results of Eq. (|6]) for different parameters. Both 
the decoherence error E T (the first term in Eq. (j6|)) and 
imperfect blockade error Eb (the second term in Eq. (|6"|)) 
agree well with the numerical results in the small gate er- 
ror limit. The total gate error with all the error sources 
is 6.5% in agreement with the atom loss corrected fi- 
delity of F = 0.92 and the simple gate error analysis 
reported in [I?) ■ As shown in Table [Til the two main er- 
ror sources limiting the gate fidelity are the spontaneous 
emission from state \p) and imperfect Rydberg excita- 
tion and blockade due to finite atomic temperature (not 
accounting for the atom loss before the CNOT pulses, 
the imperfect optical pumping and other losses that are 
independent of the CNOT pulse sequence) . 

The comparison of Monte Carlo master equation sim- 
ulations with experimental data shows that the error 
sources listed in Table [II] are able to account for mea- 
sured results with an accuracy of about 1% as regards 
the CNOT truth table and about 5% as regards the Bell 
state fidelity. The Bell state fidelity is much lower than 
that of the CNOT truth table due to dephasing that oc- 
curs while the control atom is Rydberg excited. As has 
been discussed in [ill [24| the main sources of the Ryd- 
berg dephasing are magnetic field noise and Doppler ef- 
fects. In order to significantly reduce these errors it will 
be necessary to work with colder atoms, less magnetic 
field noise, and faster Rydberg excitation pulses. 



target atom sites. Here, we treat the ground Raman pulse 
as an ideal unitary operation since it is substantially less 
sensitive to the various error sources than operations in- 
volving Rydberg states. 

The probability truth table error is defined by E = 
1 — Tr[pj^ cal p ct ], where p ct is the numerical simulation 
result. We should point out that the switching time £ gap 
is short enough that it has little effect on the CNOT 
truth table fidelity, but has a strong effect on the entan- 
glement fidelity because the motion of Rydberg excited 
atoms between excitation and dcexcitation pulses leads 
to a stochastic phase that degrades the entanglement as 
was pointed out in [l5| . Finally, we average over 100 evo- 
lutions of the master equation, and the final results are 
shown in Table [TlJ 

To model the entanglement, we prepare the control 
atom in state |c) = -^(|0) + and target atom in 

state |i) = |1). Then we follow the same approach as 
for the CNOT truth table error analysis by solving the 
time evolution of the master equation (|1 3|) . From the 
final density matrix after the pulse sequence, we can ex- 
tract the Bell state fidelity F defined as F = Tr[/9j^ eal p ct ], 
where pidoai is associated with the maximally entangled 

Bell state \Bi) = ^(|00) + |11)). We then average over merical simulations of QPT with maximum likelihood 



B. Simulated Quantum Process Tomography 

The experimental results obtained to date from Ryd- 
berg blockade experiments on pairs of atoms are far from 
predicted error thresholds for a practical fault-tolerant 
quantum computer which range from 10~ 4 — 10 -2 in dif- 
ferent models [l9l . [20| . In order to characterize more 
completely the fidelity and usefulness of Rydberg block- 
ade for quantum computing applications we need to per- 
form Quantum Process Tomography (QPT) of the 
Rydberg blockade mediated quantum blackbox process. 
QPT has been demonstrated with several different physi- 
cal systems including linear optics [3Cj| , trapped ions 31 [ , 
and superconducting circuits [32| . Here, we perform nu- 
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FIG. 4. (color online) Comparison between intrinsic gate er- 
rors from Eq. (JSJl for different parameters and numerical sim- 
ulations. In the first column (a, c, e, g): Bo/27r = 10 MHz, 
fio/27r = 1.15 MHz, ro = 300 fj,s; in the second column (b, d, 
f, h): B /2tt = 10 MHz, C1 /2tt = 0.575 MHz, r = 300 ps. 
E T is the decoherence error due to finite Rydberg lifetime (the 
first term in Eq. (|SJ) and Eb is the imperfect blockade error 
(the second term in Eq. (©). The black dotted points are the 
results of the numerical simulations and the blue lines are the 
results from Eq. © with the same parameters. 



estimation of tomographically reconstructed density ma- 
trices for the Rydberg-blockade Cz gate. We limit 
the simulations of intrinsic errors to the simpler Cz gate 
since it has been demonstrated [HJ that the additional 
single qubit pulses needed to implement a CNOT can be 
performed with errors at the ~ 10 -4 level. 

Since our goal is to determine the minimum possible 
gate error that can be reached using Rydberg blockade 
we only account for intrinsic gate errors as described in 
Sec. IIIIl and assume all additional technical errors are 
negligible. This corresponds to a situation where the 
atoms are cooled to their motional ground state so there 
is no Dopplcr dephasing during Rydberg excitation, po- 
sition dependent variations in Rabi frequencies, or AC 
Stark shifts. We assume there is no spontaneous emission 
from the intermediate \p) level during Rydberg excita- 
tion. This could be achieved using one photon excitation 



of Rydberg \p) states, or by using sufficient laser power 

|n+1,r> J. Bn . n+1 

...L B n ,n 



|n,r> 
|n-1,r> 
|n-2,r>- 



co 10 



Bnn-1 



control 



-t B n ,n-2 

target 



FIG. 5. (color online) Level diagram accounting for neighbor- 
ing Rydberg levels of control and target atoms. The Rydberg 
laser is tuned to excite |1) to \n,r). The red dashed lines 
show the energy that qubit state |0) is excited to, and B n n / 
is the positive blockade shift between a control atom in \nr) 
and target atom in \n'r). 



to detune very far from the intermediate \p) level. We 
also assume that dephasing due to time varying magnetic 
fields is negligible. 

Accounting only for intrinsic gate errors the analytical 
estimates of Sec. IIIIl show that E ~ (Br)~ 2 / 3 . At room 
temperature the Rydberg lifetime scales as r ~ n 2 with 
n the principal quantum number and in the heavy alkali 
atoms Rb and Cs the van der Waals Blockade interaction 
scales as[34| B ~ n 12 . Thus we expect the gate error to 
scale as E ~ n -28 / 3 so that choosing large n should give 
arbitrarily small errors. 

> 



This argument breaks down when n ~ 100 since the en- 
ergy spacing of levels n and nil becomes comparable to 
B or cjio, as shown in Fig. [5j This puts a limit on the ef- 
fective blockade shift that can be achieved at large n and 
limits the error floor. For the Cz pulse sequence acting 
on the four possible input states in the computational ba- 
sis {|00), 1 01) , 1 10), |11)} excitation is blocked three times 
due to wio, once due to wio + B and once due to B alone. 
It is necessary to choose the Rydberg level spacing and 
blockade shift such that the excitation suppression is as 
large as possible for all three cases. The above description 
is valid for Rydberg p 3 / 2 and d 5 / 2 states since by using 
er + polarized light the qubit states arc only coupled to 
these Rydberg states [33. For S\i2 states the situation 
is worse since the Rydberg lasers simultaneously couple 



to both 



Sl/2 



and d 



3/2 



d 5 / 2 states. 



In order to quantitatively account for coupling to more 
than one Rydberg level we have extended the basis used 
for simulations to the set {|0), \g), |1), |r 2 ) |ri), \r)}, where 
|ri)) I ^2) are additional Rydberg levels. Finding optimal 
states is now a multiparameter optimization problem. 
Details of how this is done, and the parameters of the cho- 
sen s 1 / 2 ,P3/2 an( f d 5 /2 states, are given in Appendix |A"1 In 
this extended basis, but without the \p) level, the Hamil- 
tonian and Liouvillc operators corresponding to Eqs. (|14p 
are 
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(15a) 



(15b) 



with p r T = Prr + Pnn + Pr 2 r 2 the total Rydberg excited 
population. The two-atom operators are 



H ct = H c ® I t + I c <g> H t + B ct , 



(16) 



C ct = C c <X> It + I c (8 A ; ^(c/t) are 6x6 identity 
matrices, and B ct = 7idiag[02i, B rs , r2 , B rari , B r2r , 0, 0, 0, 
B ri r 2 , B riri , B rir , 0, 0, 0, B rr2 , B rri , B rr ] is a Rydberg block- 
ade matrix where O21 is a 21 zero list. We assume that 
the Rydberg states decay directly back to the 8 ground 
sublcvels of Rb with equal branching ratios of 1/8. For 
Cs atoms the factors of 1/8, 3/4, 1/8 on the diagonal of 
(|15bp become 1/16,7/8,1/16. The Rabi frequencies for 
Rydberg excitations to states |r), |ri), \r%) are taken to 
be equal = 0' = O" for Rydberg np and nd cases. For 
Rydberg ns states we have |r) = \ns), |n) = \n — l,s), 
|r 2 ) = \n-2,d) for which O' ~ 0,0" = 1.310. Further 
details are given in Appendix|XJ For simplicity we assume 
the decay rate 7 r is the same for all Rydberg levels. It 
is straightforward to include state dependent decay rates 
in the code, but this has a negligible impact on the re- 
sults since there is very small excitation of the secondary 
Rydberg states. For simplicity we use the decay rate of 
the targeted Rydberg state for all states. 

As discussed in Ref . [36| , there is no universally agreed 
upon measure for comparing real and idealized quantum 
processes [37j . A widely used measure of quantum process 
fidelity is the trace overlap fidelity Fo, or error Eo = 
1 — Fo which are based on the trace overlap between 
ideal and experimental (in our case simulated) x process 
matrices. Another error measure Ed is defined as the 
trace distance between the ideal and simulated matrices. 
We have quantified process errors using the trace overlap 
and trace distance as 



Er 



1 - Tr 
1 



\J VXsimXidA/xI 



-Tr 



V (Xid - Xsim) f (Xid ~ Xsim) 



(17a) 



(17b) 



where Xid is the ideal process matrix and Xsim is the 
simulated physical %-matrix found from QPT accounting 



for intrinsic gate errors as described by Eqs. (|15I16[) . We 
use a maximum likelihood estimator to extract a physical 
X matrix from the QPT simulations [30j . 

In Table Hill we present the errors found from simulated 
QPT for the atomic states in Table ITVl The process to- 
mography errors tend to be 5-10 times larger than E c ^ 
which are the errors estimated in Appendix [A] for two- 
qubit product states in the computational basis. This is 
to be expected since the analytical estimates are derived 
from the probabilities of the gate succeeding, and do not 
account for output state phase errors. The trace loss 
quantifies the population in states outside the computa- 
tional basis at the end of the gate sequence. These errors 
are due to spontaneous emission from Rydberg states and 
imperfect blockade which leaves atoms Rydberg excited 
at the end of the gate. Trace loss errors account for 
about half of the process error. While the process error 
based on trace overlap Eo is less than 0.003 for all states 
listed, the error as measured by the trace distance £b is 
significantly larger. 

Some insight can be gleaned into why the trace dis- 
tance gives larger errors than the trace overlap as follows. 
The Jaksch et al. pulse sequence produces an imperfect 
Cz gate which can be written in the computational basis 
as 



C z = 



( 1 
0-10 
0-10 

V 



\ 



oi<t> 



(18) 



where <p = 7rO/2B is a small phase error in the strong 
blockade limit. As an example Fig. [6] shows the ideal 
%-matrix and the difference between the ideal and simu- 
lated x- m atrices for the standard pulse sequence leading 
to (TT8"|) . From Fig. [HJc) and (d), we can see that the error 
in the imaginary part of the x-matrix is much larger than 
that in the real part of the x-matrix. This is due to the 
fact that the real part is related to the amplitude errors 
for Rydberg blockade which arc proportional to <j> 2 , but 
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TABLE III. Gate errors from simulated QPT for several Rydberg states of Rb and Cs. The reported errors are -B c b, the 
analytical estimate found in Appendix [A] using computational basis states, trace loss, which is the sum of populations outside 
the computational basis at the end of the gate sequence, Eo(Eo) trace overlap(distance) errors from Eqs. (|17|l using the Jaksch 
et al. pulse sequence (Fig. i an d Eb the trace distance error found using the modified pulse scheme of Fig. [3b. For the s 
and p states we used the optimal Rabi frequencies found in Appendix A. For the d states we found that the Eq, E d errors were 
reduced by about 25% by using Rabi frequencies about 15% lower than those estimated in the Appendix. 





87 Rb 


Cs 


Rydberg state 




76p 3/2 


124p 3/2 


123d 5 /2 


82s 1/2 


70p 3/2 


H2p 3 /2 


H2d 5/2 


80s 1/2 


Rabi frequency £l/2ir 


(MHz) 


38.5 


16.3 


15.3 


19.2 


47.1 


19.5 


20.4 


21.4 


E c b 




.00015 


.00013 


.00016 


.00032 


.00013 
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FIG. 6. (color online) Ideal Xid matrix and the difference 
between ideal Xid and simulated Xsim matrices for QPT of 
the Cz gate using the Jaksch, et al. pulse sequence[lj for 
87 Rb 112p 1/2 states with = 27.6 MHz. (a) Real part 

of Xid, (b) Imaginary part of Xid, (c) difference of real parts, 
R- e [Xid — Xsim], and (d) difference of imaginary parts Im[xid — 

Xsim] • 



the imaginary part is related to the phase error which 
scales as 4>. 

Following the procedure for QPT in Ref. @ we can 
calculate the ideal and imperfect x- ma triccs from (fT5|) . 
and using Eqs. (|17p wc find the trace overlap and trace 



E = - [1 - COS(0)] - , 

£b = -ysin(0/2)--^~-. 



(19a) 
(19b) 



We see that Eo ~ (j> 2 while Er> ~ 4> which verifies that 
the trace overlap is not sensitive to the imaginary part 
of Xj which has linear phase errors, whereas the trace 
distance is sensitive to these errors. 

Using the modified sequence of Fig. [3]3 we can correct 
the leading order linear term in the phase error. Doing so 
has negligible effect on the trace overlap since it is only 
sensitive to amplitude errors at 0(</> 2 ). We do not report 
trace overlap errors for the modified pulse sequence in 
Table IIIII since they are unchanged. However there is a 
large reduction in the trace distance error using the mod- 
ified pulse sequence as can be seen from the values of E' D 
in the last row of the table. It has been common practice 
in experimental studies of quantum gate process fidelity 
to use the trace overlap as a reliable measure of the gate 
fidelity. The results shown in Table UTTI highlight the fact 
that the trace overlap may give an overly optimistic view 
of the gate performance, since the trace distance gives 
larger errors. Identifying what type of errors are present 
and finding ways to minimize them, as we have done here 
using a modified pulse sequence, is facilitated by checking 
several error measures. 

Finally we note that besides the intrinsic gate error 
sources, the dipolc-dipole interaction will cause a mo- 
mentum kick to both atoms [l| which can excite a trap 
state without changing the internal state of the atoms 
when they arc in Rydberg states. The perturbative tran- 
sition probability is bounded by pk < (3?7fi 2 Ai/8B) 2 /2 
with At = 1~k j£l and r\ = a/d -C 1, where a is the initial 
width of the atomic wave function determined by the trap 
and d is the trap separation. For typical gate parameters 
and rj = 1/50 wc find that p k < 3 x 10~ 7 for Rb Rydberg 
124p 3 / 2 states which is much smaller than the QPT er- 
rors in Tabic UTTI Thus errors due to momentum transfer 
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between Rydberg excited atoms have a negligible effect 
on the gate fidelity. 



V. DISCUSSION 

An important motivation for performing detailed cal- 
culations of gate errors is to determine if Rydberg block- 
ade could be used to build a fault tolerant, large scale 
quantum computing device. In order to answer that 
question it is necessary to make a connection between 
the fidelity measures, and the error limits for fault tol- 
erant architectures that have been calculated theoreti- 
cally. Threshold calculations typically proceed by assum- 
ing that an ideal unitary operator describing the time 
evolution of the quantum circuit is corrupted by an error 
operator U el with a small probability e, < e < 1 so that 
the actual gate is described probabilistically as 

U = (1 - e)U id + eU ei . 

Depending on what assumptions are made about the 
types of error operators U er that may occur, possible cor- 
relations between errors at different sites, and the over- 
all system architecture, different threshold values eth can 
be found. Provided e < eth it is in principle feasible 
to build an arbitrarily large quantum processor. Cal- 
culations that make a minimum number of assumptions 
about U er result in very low thresholds, e t h ~ 10~ 5 [38|. 
Other calculations that make more restrictive assump- 
tions result in higher thresholds. For example eth > 0.03 
in a model where the U el are Pauli operators . 

In order to relate the fidelity measures to thresholds for 
fault tolerance it is necessary to make explicit the connec- 
tion between the process fidelity and the error probability 
e, as has been done for photonic quantum gates [3!|. We 
can estimate the lower bound on e given a process fidelity 
as follows. Replace the process matrix in Eqs. (jTTj) by 
Xsim -> (1 _ e )Xid + eXcr where Xer is the process matrix 
corresponding to the operator U cr . 

To derive lower bounds on e we substitute the modified 
expression for x s i m into (|T7b ) to get 



Er 



eTr 



\J V (Xid - X^)Xid\/(Xid - Xer) 



0(e 2 



The right hand side is maximized for Xor = — Xid an d 
assuming small e we find Eo < 2e. We can therefore 
use the trace overlap to bound the error probability from 
below according to 



Er 



(20) 



Following the same steps for the trace distance gives 



e D > E L 



(21) 



These bounds result from assuming e -C 1 and a worst 
case error process with Tr[xidXcr] = — 1- Our calcu- 
lated fidelity errors given in Table Hill arc 0(1O~ 3 ) and 



we therefore have placed lower bounds on e of 0(1O -3 ) 
which is below the threshold for some fault tolerant ar- 
chitectures. 

Unfortunately this does not prove fault tolerance. We 
have bounded e from below, but the actual e for our gates 
may be higher. In addition, threshold calculations make 
assumptions about the types of errors that may occur, 
whereas our calculations of process fidelities are based on 
an independent physical model of the gate. In order to 
claim fault tolerance we would have to verify that the er- 
rors occurring in our simulations arc compatible with the 
assumptions made in the threshold calculations. This has 
only been attempted for linear optics quantum gates [39| 
and is beyond the scope of the present paper. All we 
can say based on the results obtained here is that it is 
plausible that the fidelity of Rydberg blockade gates is 
sufficient to meet the threshold for fault tolerance in an 
appropriate architecture, but this has not been explicitly 
demonstrated. 

In conclusion wc have performed a detailed analysis 
and numerical simulation of our recent demonstration of 
a Rydberg blockade mediated CNOT gate between two 
individually addressed neutral atoms. Good agreement 
between the model and experimental results allows us 
to identify the leading error sources limiting the CNOT 
truth table fidelity as imperfect state preparation, spon- 
taneous emission from the intermediate state during two- 
photon Rydberg excitation and imperfect Rydberg exci- 
tation and blockade due to variations of the atomic po- 
sition at finite temperature. The fidelity of entangled 
Bell states created so far with Rydberg blockade is pre- 
dominantly limited by ground-Rydberg dephasing due to 
Doppler broadening and magnetic field noise. 

We have also found intrinsic error limits for Rydberg 
states which are accessible by one or two photon excita- 
tion through dipole allowed transitions. We show that 
the gate error cannot be made arbitrarily small by ad- 
dressing higher lying Rydberg levels due to off-resonant 
coupling to neighboring levels which reduces the blockade 
effect. We identified the optimum blockade strength in 
the presence of neighboring Rydberg levels and showed 
using simulated QPT that for both 87 Rb and Cs atoms 
states can be found with process errors below 0.002, pro- 
vided we use a modified pulse sequence to correct small 
phase errors. The phase error correction assumes that 
the atoms are in the motional ground state of the op- 
tical traps. Our identification of optimum parameters 
including coupling to neighboring Rydberg levels is only 
approximate and it may be possible to further reduce the 
gate error with a more extensive parameter search. 

While we have focused on the Rydberg blockade 
mechanism, the direct-interaction Rydberg phase gate[l[, 
which uses simultaneous excitation of both atoms to a 
Rydberg level, may also be a route to high fidelity oper- 
ation. Recent analysis of this gate using optimal control 
theory has identified parameters for which the gate er- 
ror approaches 10 -3 although a rigorous process error 
was not calculated [i^ . We note that also the phase gate 



12 



which operates with the ordering fi 3> B will be subject 
to a limit on how high n can be due to off-resonant ex- 
citation of neighboring Rydberg levels, as illustrated for 
the blockade gate in Fig. [5] 

Our error results assume operation in a room temper- 
ature environment. The lifetimes of the P3/2 states in- 



crease by about a factor of four in a 4K He crvostat[41 



which would result in a reduction of the gate error by 
a factor of ~ 2.5 to a level below 10~ 3 . Even lower 
error levels could in principle be reached using circu- 
lar Rydberg states that have orders of magnitude longer 
lifetimes, although there are serious technical challenges 



connected with high fidelity excitation and de-excitation 
of these states. 
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Appendix A: CNOT truth table error estimates 
with multiple Rydberg levels 

In this appendix we present analytical estimates for 
the CNOT truth table including off-resonant excitation 
of multiple Rydberg levels. These estimates were used to 
find parameters for the process tomography calculations 
in Sec. El 



1. Rydberg np 3 / 2 states 

Alkali atom npi/ 2 or np 3 / 2 levels can be reached by 
one-photon excitation from the ground state. The fine 
structure splitting of high n Rydberg levels is relatively 
small, only 94 MHz for the Rb lOOp states. This small 
splitting is problematic since resonant coupling to say 
100pi/2 with a Rabi frequency of ~ 10 MHz would give 
errors at the 0.01 level due to off-resonant coupling to 
100p 3 /2- We therefore assume that the qubit state |1) is 
mapped onto the stretched ground state /+ = I + 1/2, 
rrif = f before and after Rydberg operations. When 
the stretched state is excited with a+ light angular mo- 
mentum selection rules prevent coupling to np\/2- Wc 
therefore only account for coupling to np^/ 2 ,rnj =3/2 
states. 

Referring to Fig. [7] wc assume the spacing between 
neighboring levels Lo nn _i 1 satisfies w n , n -i > uj w . This 
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FIG. 7. (color online) Off-resonant Rydberg excitation for 
p states with u> n , n -i > Wio- The filled circles indicate the 
energies where the atoms are excited to with the long-dashed 
red lines corresponding to excitation of an atom starting in 
|0>. 



corresponds to Rydberg levels with n ~ 100 for the heavy 
alkalis. The leading contributions to blockade errors for 
the computational basis states are 
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The average blockade error is Eb = \{Eqq + £m + 
Ei + En). We introduce two dimensionless parameters 
a = u) n>n -i/u)iQ and b = B„. n /c<Jio. Parameter a takes on 
discrete values as a function of Rydberg level n while b 
can be adjusted to minimize the error at fixed n by chang- 
ing the interatomic separation d. The blockade interac- 
tion between levels of different n is also a function of d, 
but to a good approximation we can put B„. n _i = 6'B„. n , 
with b' a constant independent of d. We neglect contri- 
butions from coupling to level n + 1 since the solutions 
found below have b ~ 0.5 and the average error from n+1 
states only contributes at the 10% level. 

The blockade error is minimized by ensuring that all 
undesired excitations are detuned as much as possible. 
This corresponds to w n .„_i very large but this is not a 
useful solution since it implies small n and large spon- 
taneous emission errors. In order to find a reasonable 
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FIG. 8. (color online) Blockade error scaling for p states with 
w nn _i > ujio as a function of a, b with Q/ww = 0.005 and 
b' = 1. 



value for n consider Fig. [5] which shows the blockade er- 
ror for selected values of a, b. We see that for a ~ 2.5 and 
b ~ 0.5 the error is not far from the minimum possible. 
The scaled blockade shift b could be made larger, but this 
would require very small values of the separation d which 
implies difficulty in individual addressing of the atoms. 

The above conditions are matched quite closely for 
87 Rb (Cs) using np 3 / 2 states with n = 76(70). For 
87 Rb with n = 76 at d = 1.8 /im we have 1^76,75 = 
2tt x 17.0 GHz, a = 2.49, b = 0.51, b' = 1.07 which'gives 
E B = 0.63 x 10~ 4 at VL/lo w = 0.005. For Cs with n = 70 
at d = 1.35 /im we have W7o,69 = 27rx 23.0 GHz, a = 2.50, 
b = 0.48, b' = 1.23 which gives E B = 0.66 x 10~ 4 at 
fi/wio = 0.005. 

We can estimate the CNOT truth table error aver- 
aged over the computational basis states using the same 
procedure as in Sec. IIII1 Including the spontaneous emis- 
sion errors from Eq. (|6]), neglecting corrections of order 



f2 2 /cj 2 , n 2 /B 2 , the average error is 
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where we have written Eb = EbqQ 2 - The optimum Rabi 
frequency is f2 pt = (if) 1 3 ^ Eg ^ T y/ 3 which gives the 
minimum error for the computational basis states 
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filled circles indicate the energies where the atoms are excited 
to with the long-dashed red lines corresponding to excitation 
of an atom starting in |0). 



9 



6 5 



0.3 0.4 0.5 0.6 0.7 
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For 87 Rb 76p 3/2 with r = 223 fis we find fl opt = 2tt x 
38.5 MHz (fi ptMo = 0.006) and E ch = 1.5 x 10~ 4 . For 
Cs 70p 3/2 with r = 211 /is we find rj opt = 2rr x 47.1 MHz 
(Qopt/wio = 0.005) and £; c b = 1.3 x 10" 4 . 

It is also possible to consider states with higher n such 
that uj n ,n-2 > Wio > 0J n> „-i. In this case, which is shown 
in Fig. [31 we must include off-resonant coupling to a 
third Rydberg level n — 2. The effective blockade shift 
is now smaller than for the lower n states but there is 
the advantage that the Rydberg lifetime is longer. The 
blockade errors for the computational basis states are 
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The average blockade error Eb 
is shown in Fig. [TU1 

The error is minimized for a ~ 0.57 and b ~ 0.28. 
These conditions are matched quite closely for 87 Rb (Cs) 
using np 3 / 2 states with n = 124(112). For 87 Rb with 
n = 124 at d = 4.5 /im we have W124.123 = 2tt x 3.73 GHz, 
wi24,i22 = 2tt x 7.55 GHz, a = 0.55, 6 = 0.29, V = 1.05, 
6" = 1.06 which gives E B = 3.2 x 10" 4 at fi/wi 
0.1)0"). For Cs with n = 112 at ti = 3.2 /xm we have 
Wii2,ui — 27r x 5.23 GHz, ^112,110 = 27r x 10.6 GHz, 
a =' 0.57, b = 0.29, b' = 1.19, V' = 1.18 which gives 
E B = 3.4 x 10~ 4 at fi/wio = 0.005. 

Using Eq. (JXTJ we find the following CNOT truth 
table error estimates. For 87 Rb 124p 3 / 2 with r = 616 /is 
16.3 MHz (fiopt/wio = 0.002) and 



we find f2 pt 
^cb = 1.3 x 
we find 0, o 



For Cs 112p 3 / 2 with 



r = 593 ^s 
= 0.002) and 



2tt x 
10" 4 . 

s opt = 2?r x 19.5 MHz (fi op t/wio 
E ch = 1.1 x 10~ 4 . 

We see that the truth table error estimates are slightly 
less than for the lower n situation of Fig. [7] An addi- 
tional advantage of using higher n states is that the opti- 
mal blockade shift is reached with a larger d ~ 3 — 4 (im 
which may help to minimize qubit addressing crosstalk 
in an actual implementation. For convenience we have 
summarized all parameters for the states used in Table 

m 



2. Rydberg nd 5 / 2 states 

The issue of small fine structure splitting of the p states 
discussed above, also applies to the 71.^3/2,^5/2 states. 
There are two ways of avoiding this problem. As with 
excitation of the np states we may assume two-photon 
excitation from the stretched ground state with cr + light 
which only couples to nd§/2,mj =5/2 states. Alterna- 
tively if the first leg of the excitation is made via the 
Dl transition (5s]/2 — &Pi/2 m or ® 8 i/2 — 6P1/2 m 
Cs) then only the nd^/2, rrij —3/2 states can be reached. 
Since we are setting the separation d to give the desired 
blockade strength the only difference in the gate error 
using nd 3 / 2 or nd 5 / 2 states is due to differences in the 
lifetime. The lifetimes differ by only a few percent [4l[ 
and we will therefore simply consider nd^/2 states. The 
choice of optimum states and error analysis then follows 
that in Sec. IA 11 We have summarized the parameters 
for the states used in Table HVl 



3. Rydberg nsi/ 2 states 



The ns\/2 states have no fine-structure but there is an 
additional complication since it is not possible to use an- 
gular momentum selection rules to couple to nsi/ 2 states, 
but not nd 3 / 2 or nd 5 / 2 states. We must therefore con- 
sider off-resonant coupling to additional Rydberg levels. 
The situation for uj n „_i > u)\q is shown in Fig. 1111 The 
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leading contributions to blockade errors for the compu- 
tational basis states are 
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TABLE IV. Parameters of the Rydberg states used for QPT simulations. 
Bn,n-2/B™ for np3/2,nd 5( / 2 states and b" = B^ n _i/B nn for nsi/2 states. 



The scaling parameter b" is defined as b" 
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Here a/, B' refer to frequencies and couplings between 
S1/2 and o? 3 /2 states and f2' is the Rabi frequency for 
nds/2 excitation via the Dl transition from the ground 
state so that we need only consider nd 3 / 2 states. We 
introduce dimensionless parameters a = w„ in _i/a;io, 

a' = Wn,n-l/"10) a " = "n,n-2/"10; b = B n ,„/wio, 

6 = B„ )n _i/B„ )n , 6" = B nn _ 1 /B„ jn , 6"' = B n n _ 2 /B niTl , 
and c' = Cl'/Q. Using cr_|_,cr_ polarized excitation light 
for the ground-Dl and Dl-Rydbcrg transitions d = 1.31. 
Since re — ll ^ l"n,n-2l we have neglected terms due to 



coupling to n — l,d 3 / 2 in (|A2|) . Additional checks with 
this level included give not more than 5% increase in 
the error averaged over the computational states. Since 
the computational cost of adding an additional Ryd- 
berg level in the master equation simulations is large we 
have neglected this small correction and performed mas- 
ter equation simulations with the three Rydberg states, 
|r) = |n,si/ 2 ), |n) = |n- 1, s 1/2 ), |r 2 ) = In - 2,ri 3/2 ). 

Numerical checks using Eqs. (|A2IA1[) were used to 
select the states listed in Table HVl 



